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Summary 


An iterative method for wall interference assessment and/or 
correction is presented specifically for transonic flow conditions in 
wind tunnels equipped with two-component velocity measurements on a 
single interface. The iterative method does not require modeling of the 
test article and tunnel wall boundary conditions. Analytical proof for 
the convergence and stability of the iterative method is shown in the 
subsonic flow regime. The numerical solutions are given for both two- 
dimensional and ax i symmetrical cases at transonic speeds with the 
application of global Mach number correction. 

Introduction 

The recognition of the shortcomings of the classical methods for 
determining wall interference has demanded that a variety of flow 
measurements near or on the tunnel wall be made to assess wall 
interference. For example, Arnold Engineering Development Center has 
designed an interface measurement system for the 4-foot Transonic Tunnel 
and the wall pressure measurement system is being planned for 
installation in the current 12-ft Wind Tunnel restoration project at 
NASA/Ames Research Center. It is well known that classical methods 
require a mathematical modeling of the wall boundary condition which is 
a difficult task, in particular, for a ventilated wall. The 
characteristics of a ventilated wall are not precisely known and 
generally vary with Mach number, model geometry, Reynolds number and 
wall configuration. Another shortcoming is the necessity to accurately 
represent the test article beyond the simplified representation used in 
the classical approach. 

Furthermore, the development of an adaptive wall tunnel requires 
that the boundary flow measurements be an integral part of the 
procedure. The boundary measurements are designed to determine 
adaptive-wall setting for free-air or noninterference condition. 
Techniques using only the measurements to assess the level of 
interference will complement the adaptive-wall process by determining if 
a noninterference wall setting has been achieved or by correcting for 
residual interference. Therefore, several wall interference methods 
based solely on the measured flow boundary condition have been developed 
for different kinds of measurements, such as two velocity components at 
a single interface (Refs. 1-4) or one velocity component at two 
interfaces (Ref. 5). However, all methods are limited to subsonic flow. 
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The purpose of this paper is to present a method for assessing wall 
interference for subsonic as well as transonic flow regimes based only 
on two velocity distributions measured on a single interface and no 
knowledge of the test article. This method as shown in Fig. 1 is 
divided into two tasks. The first task is the determination of the 
effective (equivalent) shape of the model from the two-measured 
variables by an iterative procedure. The second task is to predict the 
wall interference, such as a global Mach number correction, based on the 
effective shape predicted from Task 1. 

In the presentation of the paper, a proposed iterative procedure to 
accomplish Task 1 is described and convergence analytically proved for 
subsonic flow. Numerical validation is demonstrated for two-dimensional 
and axi symmetrical transonic flows using a transonic small disturbance 
equation, and the global Mach number corrections are determined for Task 2. 

Iterative Procedure of Effective Shape 


To accomplish Task 1, an iterative procedure using two measured 
distributions of velocity near the wall has been developed and is shown 
in Fig. 2. The procedure consists of two basic steps in which each 
step involves obtaining a solution for the region between a measurement 
interface (y=h) and an effective shape line (y=yi) near the model. In 
the first step, the velocity component, uo(yj), is obtained from a 
solution of the region where the boundary conditions are provided by an 
initial guess of the effective shape v°(yi), and the measured velocity 
component Ut(h). With the result, u°(yi), the second step obtains the 
solution v*(yi) in the same region using the measured boundary condition 
vt(h). These two steps are repeated until convergence has been achieved 
and the effective shape, v(yi), determined. 

Defining a merit function for the difference, v n+ ^(yi) - v n (yi), 
one can relax the effective shape between iterations by the form of 


y n+l (y 1 ) = o>u ri + 1 (y 1 ) + (1 -«) y n (y 1 ) 

with co as a relaxation factor. An updated effective shape, v(y^), is 
yielded by repeating Steps 1 and 2 with this relaxed form of velocity as 
input until the given criterion of merit function is satisfied. 

Predicted Effective Shape 


The iterative procedure for determining the effective shape is 
extremely flexible and can be applied to a variety of tunnel flow 
conditions. Depending on the flow condition, the region between the two 
surfaces could be analytically or numerically solved for the appropriate 
fluid equations, such as the Prandtl-Glauert equation, transonic small 
disturbance equation, full potential equation, or Euler equations. To 
illustrate the validity of the iterative procedure, it is applied to the 
two-dimensional, subsonic, Prandtl-Glauert equation. For this linear 
problem the solution for each iterative step can be determined 
analytically. The analytical expressions of the iterative solutions were 
derived by the Fourier transform technique as 
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The convergence and stability of these iterative procedures has been 
proved analytically in Appendix A and Appendix B, respectively. The 

results for a NACA0012 airfoil at M = 0.6 in an open-jet tunnel are 

shown in Fig. 3. It can be seen that the predicted effective shape 

compares well with an independent reference calculation from the 
original airfoil profile. 

For transonic flow, the numerical solution of the transonic small 
disturbance equation was solved for the procedure as shown in Fig. 4. 
The results for a two-dimensional 12% parabolic airfoil at M = 0.8 and 
an axi symmetrical body of revolution of a 10% parabolic arc profile at M 
= 0.975 are given in Figs. 5 and 6, respectively. For both two- 

dimensional and axi symmetrical cases, the agreement between the 
predicted effective shape and the reference calculation are very good. 
Now that effective shape has been determined, the next task is to 
determine the correction. 

Global Mach Number Correction 


The concept of global Mach number correction assumes that there 
exists an equivalent Mach number in free air where the calculated pres- 
sure distribution at a corrected Mach number will match the tunnel 
pressure distribution on the model or in a surface near the model as on 
the current procedure and simulation. This equivalent or corrected Mach 
number (Me) is computed using the effective shape determined from the 
iterative procedure. Setting the free-air Mach number equal to the tunnel 
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Mach number (Mf = Mt), free-air calculations are performed exterior to 
the surface near the model. Comparisons of the pressure distribution for 
the axisymmetric case are shown in Fig. 7. Correcting the free-air Mach 
number from Mf = 0.975 to M c = 0.925 resulted in a very good match 
between the free-air pressure distribution determined at M c and the 
tunnel pressure distribution determined from the iterative procedure as 
shown in Fig. 7. Therefore, these tunnel data at M = 0.975 correspond 
to those in free air at M = 0.925. 

Concluding Remarks 

A method of the wall interference assessment/correction for 
transonic flow conditions has been developed using a two-component 
measured flow boundary condition approach successfully. The method 
avoids the difficulty of modeling the test article and the tunnel wall 
boundary conditions. The prediction of interferences includes subsonic 
and transonic flow regimes without the limitation on the selection of 
equations. 
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Figure 2. Iterative procedure for obtaining effective flow shape at 
y = yj by using measured variables, u t (h) and v t (h). 
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Figure 4. Two-variable numerical procedure to determine effective shape. 
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Appendix A. Convergence of the Iterative Procedure 

For the convenience of mathematical manipulation, the Fourier 
Transform technique is applied to prove the convergence of the iterative 
procedure. In the transform plane, the solution for u(yi) at y = yi for 
Step 1 in Eq. (1) is of the form 


_ « . 

u [n) ( p,y ) = 

1 cosh p$(h-y { ) 


P 


_ ( n) 

v (pj^tanh pP(/i- y { ) 


(A-l) 


The boundary conditions applied to obtain the solution u(n)(yi) are 
ut(h) measured and v(n)(yi) selected. For Step 2, the solution for 
v(yi) of Eq. (2) can be expressed in the transform plane as 


v (p,h) _tn> (A-21 

v {n + U (p,y.) = - ip U ip,y )tanhp$(h-y ) 

1 coshpp(h — y ) 1 1 

where u (n) (y i ) is obtained from Eq. (A-l) of step 1 and vt(hi) is 
measured. 

Substituting u( n )(yi) from Eq. (A-l) into Eq. (A-2), one obtains a 
single expression which can be used to start the iterative procedure as 
described in the previous section 
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This single expression, Eq. (A-3), can be utilized to prove the 
convergence of the iterative procedure. For simplicity the iterative 
expression, Eq. (A-3) is written in the following function form. 
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(A-3b) 


K = tanh 2 p$ {h—yj 

Before carrying successive substitution of each iteration, we can obtain 
the converged solution directi ly from Eq. (A-3a) as v(n+l) = v( n ) = v 

; = C — (A-4) 

1 - K 
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It should be noted that the solution of Eq. (A-4) in the physical plane 
is an integral equation of the first kind which is unstable and 
difficult to solve (see Appendix B). Thus, we choose Eq. (A-3) to 
obtain the present solution in an iterative fashion. 

Let the initial guess be V(°) (yi), the n-th iteration gives 


— - 1 — K n - 

V {n) = G + GK + GK 2 + ... + GK n ~ l + v to, K n = G + v to) K n 

l-K 

Since K = tanh? p(i (h-yl) < 1, 
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< n) n 
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which converges to Eq. (A-4) as we wish to prove. 


864 



Appendix B. Stability of the Iterative Procedure 

Before examining the stability of the iterative procedure, we 
investigate the explicit solution Eq. (A-4) in the physical plane. 
Equation (A-4) in the physical plane can be expressed as an integral 
equation of the first kind 
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The effects of measurement errors of ut, vt on the effective shape will 
be evaluated by the following assumed errors 


v (x,h) 
err 


1 

- sin ax 
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(B-2) 


where a is a parameter. The measurement error will vanish as a 
approaches infinity. 


The result of effective shape due to vt prr of Eq. (B-2) obtained 
from Eq. (B-l) is 


v (x) = — sin ax- cosh a$h (B-3) 

err a 

As the vt err decreases to the vanishing small as a + «>, the error of 
effective shape is magnified to an arbitrarily large value 
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This demonstrates that the formulation of Eq. (B-l) is not stable, since 
a small measurement error in the boundary flow quantities will result in 
the effective shape of Eq. (B-4) growing to an arbitrarily large amount. 
On the other hand, we will demonstrate that the iterative procedure is 
stable as follows. 


We assume that the measurement error of vt and ut of the forms 


v (x,h) = - cos ax 

err a 
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The error of the effective shape obtained from Eqs. (1) and (2) is: 
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when the measured error vanishes as a - the error of effective shape 
Eq. (B-6) will also vanish. This indicates that the iterative procedure 
is stable as the measured input is perturbed. 
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